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The motion of a spherical dust cloud is described by the Lemaitre- Tolman-Bondi solution and is 
completely specified by initial values of distributions of the rest mass density and specific energy of 
the dust fluid. From generic initial conditions of this spherically symmetric collapse, there appears 
on : a naked singularity at the symmetric center in the course of the gravitational collapse of the dust 

On . 

cloud. So this might be a counter example to the cosmic censorship hypothesis. To investigate the 
genericity of this example, we examine the stability of the 'nakedness' of this singularity against 
odd-parity modes of non-spherical linear perturbations for the metric, i.e., linear gravitational waves. 
We find that the perturbations do not diverge but are well-behaved even in the neighborhood of the 
central naked singularity. This means that the naked singularity formation process is marginally 
■ stable against the odd-parity modes of linear gravitational waves. 
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O ■ !• INTRODUCTION 

■ The singularity theorems revealed that the occurrence of singularities is a generic property of space-time in general 
relativity However, those theorems say nothing about the detailed features of the singularities themselves; 

— «. i for example, we do not get information from those theorems about whether the predicted singularity is naked or 
not. Naked means that the singularity is observable. A singularity is a boundary of space-time. Hence, in order to 
obtain a solution of hyperbolic field equations for matter, gauge fields and space-time itself in the causal future of a 
naked singularity, we need to impose a boundary condition on it. However, we do not yet know physically reasonable 
boundary conditions on singularities and hence to avoid this difficulty, the cosmic censorship hypotheses (CCH) 
proposed by Penrose UJ^] are often adopted in the analysis of the physical phenomena of the strong gravitational 
k>( | fields. 

There are weak and strong versions of the CCH. The weak CCH states that a singularity is covered by an event 
horizon and never observed by anyone included in the causal past of future null infinity (not globally naked) while the 
strong CCH says that nobody can observe a singularity (not locally naked) . However, the validity of the CCH is one 
of the most important open questions in classical general relativity. No one has ever proved that these hypotheses 
hold. On the contrary, some researchers found, analytically or numerically, that there are solutions of the Einstein 
equations which have naked singularities. If these naked singularities are physically realizable, then we could be in an 
embarrassing situation because an important assumption in theorems on the nature of a black hole is violated. In the 
vicinity of a singularity, quantum effects on the space-time will play an important role and therefore if the existence 
of naked singularities is generic, in order to understand the nature of a black hole, we might need the knowledge of 
the quantum theory of gravity even if the black hole is a classical entity. 

In the last two decades several researchers have shown that the Lemaitre- Tolman-Bondi (LTB) space-time [[||7]] 
is one of the candidates for a counter example to both versions of the CCH. This space-time describes the motion 
of a spherically symmetric inhomogeneous dust cloud and is completely specified by initial values of the rest mass 
density and specific energy of the dust fluid. Eardley and Smarr showed that the central singularity of the LTB 
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space-time can be shell-focusing and naked in the case of marginally bound collapse M. Christodoulou showed that 
the same is true also for the bound case ||. Newman generalized Christodoulou's analysis to cover a larger class of 
LTB space-times 119]. Joshi and Dwivedi carried out studies of a much more general class of solutions in which a 
conical singularity (but not curvature one) was allowed in the initial configuration and showed that the formation 
of a central naked singularity is a general feature for a very wide range of initial data in the LTB space-time pi) . 
These results are summarized as follows; in this space-time, a naked singularity appears from generic initial data for 
spherically symmetric configurations of the rest mass density and specific energy of the dust fluid. 

In order to recognize this example as a serious counter example to the CCH, we should examine its genericity. That 
is to say, there is a possibility that the naked singularity is due to physically unrealistic conditions, e.g., assumptions of 
the spherical symmetry, dust matter and so on. Shapiro and Teukolsky studied evolution of collisionless gas spheroids 
by fully general relativistic simulations [12| . They found that prolate spheroids with sufficiently elongated initial 
configurations and even with some angular momentum, may form naked singularities. Ori and Piran numerically 
examined the structure of self-similar spherical collapse solutions for a perfect fluid with a barotropic equation of 
state They showed that there is a globally naked singularity in a significant part of the space of self-similar 

solutions. Joshi and Dwivedi analytically investigated the self-similar spherically symmetric collapse of a perfect 
fluid with a similar equation of state [ [l4| and further the naked singularity produced by the gravitational collapse of 
radiation shells Jl5| and of more general matter [H. 

In this article, we concentrate our attention on the issue whether the spherical symmetry is essential to the occurrence 
of the shell focusing naked singularity. For this purpose, we consider odd-parity modes of non-spherically symmetric 
perturbations in the marginally bound LTB space-time and examine the stability of the 'nakedness' of that naked 
singularity against those linear perturbations. As for the non-spherically symmetric collapse case, Joshi and Krolack 
revealed that a naked singularity appears also in the Szekeres space-time with an irrotational dust matter |l7[ . Since 
the odd-parity perturbations correspond to the rotational motions of the dust fluid and of the space-time itself, our 
analysis will give a new insight into the formations of naked singularities in non-spherically symmetric space-time. 

To decouple physical effects from gauge or coordinate ones, we adopt the gauge-invariant formalism formulated 
by Gerlach and Sengupta [fLsll for general spherically symmetric space-times. Here we consider only the metric 
perturbations, i.e., linear gravitational waves. Using this formalism, we obtain a single decoupled partial differential 
equation for a gauge-invariant variable corresponding to the odd-parity metric perturbations of the LTB space-time. 
We analyze this equation numerically by use of single null coordinates, which was adopted by Goldwirth and Piran Jl9| ] 
for the numerical study of spherical collapse of a massless scalar field. Then we shall discuss the stability of the LTB 
space-time with a central naked singularity from the results of these analyses. A naked singularity is interpreted to be 
unstable, if perturbations tend to diverge as they approach the naked singularity and the Cauchy horizon associated 
with it. If such a behavior is found, it means that the perturbations will destroy the Cauchy horizon and change the 
causal structure of this space-time. Waugh and Lake examined the stability of the central naked singularity of the 
LTB space-time against perturbations of a massless field by the use of the high-frequency (eikonal) approximation 
and with the assumption of no back reaction of the massless field to the space-time geometry p0|] . Their analysis 
revealed that the formation of the central naked singularity is stable within the validity of their approximation. In 
contrast to the analysis by Waugh and Lake, however, the effect of the finite wavelength of the perturbations and 
the non spherically symmetric dynamics of the space-time itself are taken into account up to the linear order in our 
present analysis. 

This paper is organized as follows. In Section 0, we briefly describe the LTB space-time. In Section [II, we derive 



the basic equations for perturbations in the LTB space-time and then give regularity conditions for the perturbations 
at the symmetric center, which are significant for our stability analysis of the central naked singularity formation. 
In Section IV, we present the expressions for the perturbations of the Riemann tensor. We show the numerical 
procedure and results for the marginally bound LTB space-time in Section Section VI is devoted to the summary 



and discussions. We adopt the geometrized units, c = G = 1. The signature of the metric tensor and sign convention 
of the Riemann tensor follow Ref . ^lj . 

II. LEMAITRE-TOLMAN-BONDI SPACE-TIME 

The inhomogeneous spherically symmetric dust collapse is described by the LTB space-time. Using the synchronous 
coordinate system, the line element of this space-time is expressed in the form 

ds 2 = g^dx^dx" = -dt 2 + A 2 (t, r)dr 2 + R 2 (t, r){d9 2 + sin 2 6d<j) 2 ). (2.1) 

The energy-momentum tensor for the dust fluid is 

= p{t,r)u^u v , (2.2) 
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where p(t, r) is the rest mass density and u' 1 is the 4- velocity of the dust fluid. In the synchronous coordinate system, 
the unit normal vector field to the spacelike hypersurfaces is geodesic and there is a freedom of which timelike geodesic 
field is adopted as the hypersurface unit normal. Using this freedom, we can always set = Sq since the 4-velocity 
of the spherically symmetric dust fluid is tangent to an irrotational timclike geodesic field. 

Then the Einstein equations and the equation of motion for the dust fluid reduce to the following simple equations, 

A= drR , (2.3) 

/ x 1 1 dF i r ) 

^ r) = l*iwi*' (2 - 4) 

(d t R) 2 - ffl = /(r ), (2.5) 



where f(r) and F(r) are arbitrary functions of the radial coordinate, r. From Eq.(2.4), F(r) is related to the Misner- 
Sharp mass function p^ |, m(r), of the dust cloud in the manner 

m{r)=4:TT p(t,r)R 2 dR = 4?r / p(t,r)R 2 d r Rdr = (2.6) 

Jo Jo 2 

Hence Eq.( [2.5| ) might be regarded as the energy equation per unit mass. This means that the other arbitrary function, 
f(r), is recognized as the specific energy of the dust fluid. The motion of the dust cloud is completely specified by 
the function, F(r), (or cquivalcntly, the initial distribution of the rest mass density, p) and the specific energy, f(r). 
When we restrict our calculation to the case that the symmetric center, r = 0, is initially regular, the central shell 
focusing singularity is naked if and only if d 2 p\ r= o < is initially satisfied for the marginally bound collapse, f(r) = 
p3| , p4f . For the collapse that is not marginally bound, there exists a similar condition as an inequality for a value 
depending on the functional forms of F(r) and f(r) . 

III. PERTURBATION METHOD 

The perturbation method used in this article is the gauge-invariant one which has been formulated by Gerlach and 
Sengupta [jig] for a general spherically symmetric background space-time. First we briefly review their formalism only 
for the so-called odd-parity modes. Thereafter we will apply this formalism to the LTB space-time in order to derive 
the basic equations for our analysis. 

The perturbed metric tensor is expressed in the form 

g^u = Sfiu + h^, (3.1) 

where is the metric tensor of the spherically symmetric background space-time and is a perturbation. The 
energy-momentum tensor is written in the form 

Tfo, = + 5T^ IV , (3.2) 

where T M „ is a background quantity and 5T^ V is a perturbation. By virtue of the spherical symmetry of the background 
space-time, is expressed in the form 

T^dx^dx" = f ab dx a dx b + ^f B B R 2 (t, r)dtt 2 , (3.3) 

where the sub- and superscripts, a, b, . . . represent t and r while A,B,... represent 9 and <p. The odd-parity pertur- 
bations of and ST^ are expressed in the form 

h a (t,r)^r B 

V= I fti(*,r)*p fl | , (3.4) 



Sym h 2 (t,r)xi 



AB 



t Q (t,r)^ B 

5T^= \ ii(t,r)0>r s | , (3.5) 



Sym t 2 (t,r)x" 



AB 
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where and xTab are odd-parity vector and tensor harmonics associated with the spherical symmetry of the 

background space-time [{25)] . 

We introduce gauge-invariant variables defined by Gerlach and Sengupta. The metric variables are given by 



k a = K l -R 2 d a 
The matter variables are given by the following combinations, 



L = h- -T B B h 2 . 



Then the Einstein equations lead to the equations for the metric variables as 

fc% = 16rrL (l>2), 



(3.6) 

(3.7) 
(3.8) 

(3.9) 



(R 4 W ab ) , 6 + (J - 1) {I + 2) k a = lfaR 2 L a 



(*>1), 



where W a b is defined as 



W n 



hb 
R? 



R 2 



h 
R 2 



R? 



(3.10) 



(3.11) 



\h 



and the vertical bar refers to the covariant derivative within the 2-dimensional sub-space-time (t, r). From the equation 
of motion for the matter, we get 



(i? 2 L a ) =(/-l)(/ + 2)L (Z>1). 



(3.12) 



Now we apply the above formalism to the case of the background LTB space-time. From Eqs.(2.2) and (p.q), we 
find that there is no density perturbation and that only the perturbation of 4- velocity Su^, exists; 



5u^ = (0,0,U(t,r)<P] n B ). 
Therefore the odd-parity gauge-invariant matter variables become 

Lo = pU and L\ = L = 0. 



From Eqs.(3.9) and (3.10), we obtain the equations of motion for the metric variables, 



d t (Ako) - d r (^j = 0, 
d r (R 4 ip s ) + A (I - 1) (I + 2) fc = l6nAR 2 L Q 
d t (i?V s ) + ~ (I ~ 1) (I + 2) fci = 0, 
where we have introduced another gauge- invariant variable, -0 S , defined as 



^=1 



fry 
R 2 



Eq.( p.l2[ ) becomes 

d t (AR 2 L Q ) = 0. 
This equation is easily integrated and we obtain 

dJ(r) 



AR 2 L a = 



dr 



(3.13) 
(3.14) 

(3.15) 
(3.16) 
(3.17) 

(3.18) 
(3.19) 
(3.20) 
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where J(r) is an arbitrary function depending only on r. From Eqs.(3.16), (3.17), and ( 3.2d| ), we obtain a single 
decoupled wave equation in the form 



d t (Jgdt (iZ 4 ^.)) ~ Q r (j& 9 r ( i?4 ^)) + (i " 1) (i + 2) 



-167T<9 r 



( \ dJ 



(3.21) 



The variable, ^i> s , differs from the gauge- invariant variable used in Ref. [£6| by a factor of 1/R 2 . The reason why 
we adopt tp s as a gauge-invariant variable instead is due to regularity conditions at r = 0. Further, as will be shown 
later, ip s is closely connected with the curvature tensor near the center. 

Let us consider the regularity conditions for the background metric functions and gauge-invariant perturbations at 
r = 0. Hereafter we restrict ourselves to the axisymmetric case, i.e., m = 0. Note that this restriction does not lose 
generality of our analysis. Further we consider only the case in which the space-time is regular before the occurrence 
of the singularity. This means that, before the naked singularity formation, the metric functions, R(t,r) and A{t, r), 
behave near the center in the manner 



R 
A 



R c (t)r + 0(r 3 ), 
R c (t) + 0(r 2 ). 



(3.22) 
(3.23) 



To investigate the regularity conditions of the gauge- invariant variables, k a and Lq, we follow Bardeen and Piran p7| 
The results are given by 



La- 


- L c (t)r l+1 


ke — 


- k 0c (t)r l+1 


fci - 


- k lc (t)r l + 2 



+ 0(r l+3 ), 
+ 0{r l+3 ), 
+ 0{r l+i ). 



From Eqs.( |3.18D , fl3.22| ), fl3.23D , fl3.25D and fl3.26| ), we find that ij; a behaves near the center as 



4> s 



i, sc {t)r 1 - 2 +0(r l ) 
V»-c(t)r + 0{r 3 ) 



for I > 2, 
for I = 1. 



In the case of I > 2, the coefficient, ip sc (t), is related to R c (t) and ko c (t) in the manner 

k 0c (t) 



(3.24) 
(3.25) 
(3.26) 



(3.27) 
(3.28) 



(3.29) 



From the above equations, we note that only the quadrupole mode, I = 2, of tp s does not vanish at the center. 



IV. PERTURBATION OF RIEMANN TENSOR 



In this section, we consider the perturbation of the Riemann tensor, i? M „ (T A , of the LTB space-time to investigate 
the relation between the singularity formation and the perturbations. The Riemann tensor is decomposed into the 
Ricci tensor, R^, and the Weyl tensor, 



We shall give them in the form of the components of the following tetrad basis, 

hoPte 



(*) 



'(r) 



'(0) 



1,0,0, 



1 



R 2 sin 9 
hxPr 



AR 2 sm0J , 



1 



0,0,- 
0,0,0, 



R' 2i? 3 sin 2 i 
1 



Rsm9 



(4.1) 

(4.2) 
(4.3) 
(4.4) 
(4.5) 
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where Pi(cos9) is the Legendre polynomial and the comma followed by 6 denotes a derivative with respect to 0. The 
Weyl tensor is then decomposed into the so-called electric part, E a p, and magnetic part, B a p, which are defined as 



E a [3 = Ca/j/3v6/j)6(t)) 



(4.6) 
(4.7) 



where e Mt , a /3 is the 4-dimensional skew tensor. In the background LTB space-time, the Ricci tensor has a non-zero 
value in the region of non-vanishing rest mass density, p ^ 0, through the Einstein equations and also the electric 
part has a non-zero value. On the other hand, the magnetic part is identically equal to zero in the background LTB 
space-time. However, when axisymmetric odd-parity metric perturbations exist, the Riemann tensor is perturbed and 
the magnetic part may also have a non-vanishing value. 

The perturbation of the Ricci tensor is expressed by the matter perturbation through the Einstein equations as 



8tt 



LqPl 



8tt dJ 
AR 3 ~d~r 



P 



1.0, 



(4.8) 



and the other components vanish, where we have used Eq.( |3.2d ) in the last equality. The perturbations of the tetrad 
components of the electric part arc given in the form 



S ( E (8)(<t>)) 



AR 3 



(l-l)(l + 2)k 1 +R(d t R)ips 



sin dP 



2AR 2 



(cU)fco 



(smOPifift - cos 0Pi,e) ■ 



(4.9) 
(4.10) 



and the other components vanish. The perturbations of the tetrad components of the magnetic part are obtained in 
the form 



^l(l + i)M, 



1 



AAR 3 
1 



[R 2 d r (R 2 ip s ) - A(l - + 2)ko] Pi,e, 



5(B WW ) = - 



AR 2 
1 

2AR 2 
1 

AR 2 
1 

2AR 2 



Rd r 



dtR_dtA 
R A 



fci + -AR 2 tp, 



no, I § 

Rdr[k i 



d t R d t A\ 

- - L r )ki+ AR 2 ip s 



R 



A 



d t R d t A 



1 



dtR d t A 
R 



-j- ] fci + AR £ i>, 



l(l + l)Pu 

cotep h e 
i(i + i)Pi, 



(4.11) 
(4.12) 

(4.13) 
(4.14) 



and the other components vanish. 

Now we will investigate the behavior of th e Ricc i and Weyl tensors near the center where the naked singularity 
appears. From the regularity conditions ( 3.22 )-( 3.29 ) , we can see that the perturbations of the Ricci and Weyl tensors 
obtained in the above behave near the center in the manner 



for the Ricci tensor, and 



1 

2Rj 
1 

2~m 



(i -i) 



(l + 2)k lc - R c ^k 0c 
at 



dR 

(l + 2)fci c - Rc-^-k 0c 



{Pi,e,e - cot 6Pij)r 



1 



.^(l- 1)1(1 + l)k 0c P ir 1 - 2 , 



(4.15) 

(4.16) 
(4.17) 
(4.18) 
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8(B {T . m ) — » - — {I - 1)(Z + l)/c 0c P/, 9 r z - 2 , (4.19) 

— > (Z + l)fcoc(fl,fl,9 + lPi)r 1 - 2 , (4.20) 

S(B WW ) ^-^{l + l)k 0c {cot6Pi,B+lPi)r 1 - 2 , (4.21) 

for the Weyl tensor of I > 2. For the I = 1 mode, we find 

*(i5 (r)W )— ►i^ TC r a sin0, (4.22) 

<5(-B(r)(r)) — >-4> sc rcos6, (4.23) 

$(•%)(<?)) — >• -V> sc rsin0, (4.24) 

<K-B(0)(<»)) — >-^ sc rcos6, (4.25) 

= (4.26) 



From the above equations, we see that the perturbations of the tetrad components of the Ricci and Weyl tensors, 
except for the quadrupole mode, I = 2, of the magnetic part, B a p, identically vanish at the center. This means 
that the central naked singularity formation is affected only by the quadrupole mode up to linear order. Therefore, 
hereafter we shall consider the quadrupole mode only. On the other hand, since the solution for the dipole mode, 
I = 1, is obtained analytically, we present it in Appendix although the dipole mode vanishes at the center from 
the regularity condition and does not influence the formation of the central naked singularity up to linear order. 



V. NUMERICAL RESULTS AND DISCUSSIONS 



Here we will perform numerical integration of Eq.( 3.2l] ) for the quadrupole mode, 1 — 2. As mentioned above, the 
gauge- invariant matter perturbation variable, Lq, and the perturbation of the Ricci tensor vanish at the regular c enter . 
Here we restrict our investigation to the case of no matter perturbations, namely the right hand side of Eq.(3.21) 
vanishes. The non-vanishing matter perturbation case should be investigated in future. Further, for the simplicity of 
calculations, we consider only the marginally bound cas e, f(r) = 0, as background space-time. 

By virtue of f(r) — 0, we can easily integrate Eq.(|2.5|) and obtain 



/?(/.,) = [—) [t (r) 



~i] 2/3 , 



(5.1) 



where to(r) is an arbitrary function of r. Using the freedom for the scaling of r, we choose i?(0, r) = r. This scaling 
of r corresponds to the following choice of to(r), 



to(r) 



r.3/2 



3%/F 



Here note that, from Eq.(2.3 ), th e background metric variable, A, is equal to d r R. 
Then, the wave equation^ 3.21) becomes as follows, 



1 d 2 ^ s 



1 



dt 2 {d r R) 2 Or 2 



d r R d 2 R\ dt 



{d r R) 2 V R d r R J dr 

d t d r R \ d t R fdtR^ 2 
d r R J ~R + \~R 



-4 



R 



(5.2) 



d t R d t d r R\ dtp 



R 



d r R J dt 



(5.3) 



We solve this partially differential equation numerically. In the rest of this section, we explain the details of the 
background space-time considered here, numerical methods and boundary conditions. Further we show the numerical 
results. 
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A. Background Density Profile 



The background space-time and the motion of the marginally bound dust cloud are completely determined by the 
initial rest mass density profile, p(0,r). Further as already mentioned, it should satisfy the condition d^p\ r= o < in 
order that the central naked singularity is formed. We, therefore, adopt the following initial rest mass density profile 
so that the central naked singularity appears; 

,-(0,0 = ( ^ - 2 ^ r f + {r ' r ^ ** < r < * (5.4) 
rv ' ' [ for r > rb, v ' 

where p$ is a positive constant and r b denotes the radial coordinate at the surface of the dust cloud. The total 
(gravitational) mass of the dust cloud is 

32tt 

M = m(r b ) = j^rPorl (5.5) 
The time of the central naked singularity formation is 

f = fo(0) = — L=. (5.6) 

Whether the naked singularity is global or local is determined by a non-dimensional constant po r t- It is known that 
the singularity is globally naked for sufficiently small p^r\ P,P3|. However, the critical value of p^r\ can not be 
obtained explicitly. Hence, after por% is given, we have to investigate whether the central naked singularity is global 
or local, by numerically solving the future directed null ray from the central naked singularity. Here we consider two 
cases. One is that of pqt\ = 3 x 1CP 2 , which corresponds to a globally naked singularity, and the other is that of 
Po r b = 3 x 10 _1 , which corresponds to a locally naked one. 

In the globally naked case, the initial radius of the dust cloud and the time of the central naked singularity formation 
are given by 



R(0,r b ) r b 105 



M M 32np rl 
t (Q) 105 



M 32\/67r 3 / 2 (po^) 3 / 2 
On the other hand, in the locally naked case, they are given by 

R(0,r b ) r b 



34.8, (5.7) 
= 46.3. (5.8) 



M M 
*o(0) 



3.48, (5.9) 



M ~ 1M - (5 ' 10) 



B. Numerical Procedure 



Next, we describe the procedure of our numerical calculation. We have a disadvantage when we use the (t, r) 
coordinate system, because of the restriction on the region in which we can numerically construct the solution of the 



wave equation, (3.21). Therefore, instead of the (t, r) coordinate system, we introduce a single-null coordinate system, 
(u, r'), where u is an out-going null coordinate and chosen so that it agrees with t at the symmetric center and we 
choose r' — r. We perform the numerical integration along two characteristic directions. The transformation matrix 
is formally expressed in the form 

dr' = dr, (5.11) 
du = {d t u) r dt + (d r u) t dr. (5.12) 

Because u is the out-going null coordinate, the following relation holds, 

{d r u\ d r R { ' ' 
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Using these relations, we obtain the line element of the 2-dimensional sub-space-time, (t, r), in the following new form 



^ s (2) = —o^du 2 — 2a(d r R)dudr' , 



(5.14) 



where we have introduced 



(d t u) r ' 



Furthermore, from Eqs.(pll) and ( fU2| ), we obtain 



9 r' = -S^ 9 * + d r = ( d rR)9t + d r , 

(o t u) r 



(5.15) 



(5.16) 



where we have used Eq.(5.13) in the last equality. The above equation describes that d r > is parallel to the future 
directed out-going null direction. 



By using this new coordinate system, (u, r'), Eq.(5.3) is expressed in the form 



du 



a 
R 
a 
2 



3d r R + -R(d t R)d t d r R - {d t R) 2 d r R + ^R(d r R)d?R 
J^L _ 1 (i _ dtR) 



1 

R ( 



{d r R) 
d r R 



R 



where the ordinary derivative in the left hand side of Eq.(5.17) is given by 



du 



= d„ 



dr' 
du 



d r ' = d u 



2d r R 



Or', 



(5.17) 
(5.18) 

(5.19) 



and we have introduced a new variable, <p s , given by Eq.(5.18). 

The procedure of the numerical integration is as follows. At the first step, we prepare initial data corresponding to 
imploding waves for <f> s at each grid point on the initial null hypersurface labele d by u — uq ^constant. Then, using 
this cf> s , "0s is obtained at each grid point on u — uq by the in tegra tion of Eq.( 5.18 ). At the next step, in order to 
obtain <f> s at each grid point on u = uq + Au, we integrate Eq.(5.17) by using values of 4> s and ip s on u = uq. Then, 
tp s on u — uq + Au is obtained from Eq.( 5.15 ) by using (j) s on u = uq + Am. We repeat this procedure and obtain a 
solution outside the Cauchy horizon associated with the central naked singularity. 

In the above procedure, we should impose a boundary condition on ip s at the center to perform the numerical 
integration of Eq. Q5.18 ). From Eqs.(3.27) and ( 3.29] ), in the case oil = 2 mode, ip s behaves in the manner 



+ 0(r 2 ) for 



0. 



(5.20) 



Hence, we have to numerically make ips near the center so that Eq.( [5.20" ) is guaranteed on a surface of t ^constant, 
and this leads to the boundary condition for ip s at the center. 

Here we comment on our numerical code. We compared the numerical results for the Minkowski space-time with 
the analytical solutions which will be described in the next subsection. In this case, the result produced by our code 
agreed with the analytical solution very closely. Another check we performed was to compare the numerical results 
for several mesh sizes with each other. This test confirmed that our numerical results were almost independent of the 
mesh size. 



C. Initial Conditions and Numerical Results 



The initial conditions which we consider are a Gaussian-shaped wave packet with respect to the coordinate, r' 

(r' - r' c f 



ipl exp 



2cr 2 



(5.21) 
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where tp l s , a, and r' c are constants and characterize the amplitude, width and initial position of the initial wave packet, 
respectively. The initial null hypersurface, u = uq, is chosen so that it includes a world point (t,r) = (0,0), except 
for the analysis of the scattered waves which will be discussed in this section. 

We investigate models with three different initial positions of the wave packet, i.e., r' c in Eq.( |5.2l[ ), on the initial 
null hypersurface. In Case 1, the wave packet reaches the center of the dust cloud before the formation of the central 
naked singularity. In Case 2, a significant portion of the wave packet hits the central naked singularity. In Case 3, the 
packet does not hit the central naked singularity but reaches the Cauchy horizon associated with it. Fig.[j] shows these 
situations schematically. In each case, the value of ip s at the center is plotted as a function of the coordinate time, t, 
in Fig.|| for the globally naked case and in Fig ^3] for the locally naked case. Note that it is impossible to perform the 
numerical calculation in the causal future of the central naked singularity. Therefore we plot ip s at the center only 
before the occurrence of the central naked singularity. Although such a difficulty exists, we find that violent growth 
of the amplitude of ijj s is not observed near the central naked singularity and Cauchy horizon associated with it. 

Next we show the dependence of ips at the center on the width of the initial wave packet. Fig.^ depicts ip s at the 
center for various widths of packets in Case 2. It is found that the amplitude of ip s with smaller initial width becomes 
larger at the center. The relation between the width, er, and the maximal value of |^ s | at the center is shown in FigJ^. 
We find that there is the following power-law relation 

|^|r=0 oc cr^ 3 . (5.22) 

We also observe the time dependence of ip s along the line of a constant circumferential radius outside the dust 
cloud. Since we would like to see the effect of the central naked singularity on ip s , we consider the globally naked case 
only. We set up an initial wave packet of a = 0.05rb at R — 100M on the initial null hypersurface which does not 
include the space-time point (t, r) = (0, 0) but is chosen so that the wave packet will reach the neighborhood of the 
central naked singularity. 

The results are shown in Fig.|| in which ip s at R = lOOAf is plotted as a function of t. Note that the point, 
R = 100M, is located in the vacuum region which is the Schwarzschild space-time by Birkhoff's theorem. Hence the 
value of t along the curve of R = 100M agrees with that of the usual static time coordinate of the Schwarzschild 
space-time. 

In Fig.^(a), the solid line corresponds to Case 1 while the broken line is for Case 2. The dotted line denotes the 
result for Case 3. The left-hand peaks in Fig.^(a) correspond to the initial incident waves. On the other hand, the 
right-hand peaks of Cases 1 and 2 in this figure correspond to the scattered outgoing waves. In Case 3, the right-hand 
peak does not exist and this is because, in this case, almost all portions of the incident waves enter into the Cauchy 
horizon associated with the central naked singularity and hence it is impossible to follow numerically the scattered 
waves in the causal future of the central naked singularity. Fig.^(b) shows detailed behavior of the scattered ip s for 
Case 2. It is a most important fact seen in these figures that the amplitude of the scattered waves is almost the same 
as that of the initial incident waves in Cases 1 and 2. 

In order to investigate the effect of the wavelength of ip s , we perform the numerical integration for Case 2 but with 
different initial widths of wave packets. The results of narrower (a = 0.027-f,) and broader (er = 0.25rb and 0.5r&) 
widths than the case plotted in Fig.^|(a) and (b) are shown in Fig.g(c). The narrower wave is similar to a = 0.05rb 
while the broader packets have different forms of scattered waves from the narrower one. However, in both cases, the 
amplitude of the scattered wave is not so different from the incident one. 



D. Minkowski Case 



Here we investigate the behavior of ip s in the Minkowski space-time and compare it with the results of the LTB 
space-time obtained in the above in order to reveal the effects of the existence of the dust cloud and the central naked 
singularity on the dynamics of ip s . In the Minkowski case, since R(t,r) — r, Eq.(5.3) becomes 

dfy s -d 2 r ip s = -d r *p s . (5.23) 
r 

The solution of this equation which is regular at r — is obtained in the form 

J(t-r)-f(t + r) f M(t-r) + fM(t + r) fW(t - r) - fW(t + r) 
^ = 3 ~ 5 + 3 - 4 + ~ 3 , (5.24) 

where fix) is an arbitrary function and f^ n \x) denotes the n-th order derivative of fix) with respect to x. We set 
the following initial wave packet on the null hypersurface, t = r, 
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4> s 



exp 



(r - r c ) 
2<r 2 



21 



(5.25) 



Using the above solution, we compare the evolution of wave forms in the Minkow ski s pace-time with that in the LTB 
space-time. The initial wave packet in the LTB space-time has been given by Eq. ( 5.21 ) as a function of the coordinate 
radius, r' . However, note that r' does not agree with the circumferential radius, R, in this case but in the Minkowski 
case, the coordinate radius, r, agrees with the circumferential radius, R. Since the circumferential radius, i?, is tightly 
connected with the behavior of the amplitude of the wave, we should set the sam e init ial data with respect to R both 
for the LTB and Minkowski cases. Hence first we plot the initial wave packet (5.21) as a function of R/M on the 
initial null hypcrsurfacc and then the values of a and r c in Eq.(5.25) are adjusted so that the initial wave form fits 
well with that of the LTB case. 



First we consider the evolution of ip s at the center. Using Eqs.(5.24) and (5.25), we obtain ip s at the center in the 
form 



<M',o) = 



1-Xd-rM- 



1 

4^2 



1 

4^4 



20a 4 



60cr e 



t 

r, 

2 



exp 



1 

2^2 



(5.26) 



The parameters, a and r c , in Eq.( 5.25 ) are chosen so that the initial wave packets fit well with those of the Cases 1 
and 2 of Fig^. The results are given in Fig.^(a) and (b), respectively, and in this figure, we also plot the results for 
the corresponding cases of the LTB space-time. It should be noted that there is scarcely any difference between the 
wave forms of the Minkowski and LTB cases. 

Next we consider the behavior of ip s at a finite circumferential radius which agrees with the numerical value of 
R = 100M in the LTB case. Here the wave form is obtained numerically by the same procedure as in the LTB case. 
The result is shown in Fig.||. We also plot the corresponding case of the LTB space-time in the same figure. We find 
that there is a little difference of the phase between the Minkowski and LTB cases. However, the behavior of ips in 
the Minkowski case is basically the same as that in the LTB case. The effect due to the dust cloud and the existence 
of the central naked singularity on the propagation of ip s is rather small. 

We consider the relation between the maximum value of \ip s \ observed at the cent er and the width, tr, of an initial 
wave packet in the Minkowski space-time. This relation is obtained from Eq.(5.26). The results are also shown in 
Fig.|. The power-law relation Eq.( |5.22| ) is also valid in the Minkowski case. From Eq.( |5.24| ), ip s is approximately 



proportion to 1/r 3 except for the region of r < a around the center. If the initial amplitude of the wave packet has a 
value ipi at r — r c , then the value of ip s at r = a is roughly estimated as ipi x (er/r c )~ 3 . This will be the reason why 
the relation ( 5.22| ) holds in the Minkowski space-time. As we have discussed above, ip s behaves outside the Cauchy 
horizon of the LTB space-time in ap proxim ately the same manner as in the Minkowski space-time. Therefore it would 
be also the reason why the relation (5.22) holds in the LTB space-time. 

As a result, we conclude that even in the neighborhood of the central naked singularity and of the Cauchy horizon 
associated with it, the metric perturbation, ip s , does not show any peculiar behavior. However, we should note that 
ip s does not vanish in the neighborhood of the central naked singularity although it is well-behaved. Therefore, the 
formation process of the central naked singularity is marginally stable against the odd-parity metric perturbations. 



VI. CONCLUSIONS 



We have investigated the behavior of the odd-parity linear perturbations in the LTB space-time. We have derived 
the wave equation for the gauge- invariant variable, ip s . From the analysis of the regularity for ip s and the perturbations 
of the Ricmann tensor, only the quadrupole mode, / = 2, of ip s and of the magnetic part of the Weyl tensor does not 
vanish at the symmetric center of the background LTB space-time, where a naked singularity appears in the course of 
the gravitational collapse of the dust cloud. Therefore this quadrupole mode is the most important for the stability 
analysis of naked singularity formation in the LTB space-time. Then we have performed numerical experiments on 
how a Gaussian-shaped incident wave packet behaves under this wave equation for the I — 2 mode without matter 
perturbations. From those numerical experiments, we have obtained the following results. When this wave packet 
approaches the center, its amplitude becomes larger but finite. The amplitude at the center depends on the width of 
the initial wave packet according to a power law. On the other hand, when the incident wave packet initially located 
outside the dust cloud returns back to the same circumferential radius as the initial one, the amplitude of the returned 
wave is almost equal to that of the incident one. 
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In order to reveal the characteristic effects of the LTB space-time on the behavior of ip s , we have also investigated 
%j) a in the Minkowski space-time. Then we have found that in the outside of the Cauchy horizon associated with the 
central naked singularity, the behavior of ip s in the LTB space-time seems to be not so different from that in the 
Minkowski space-time at least except for the extreme neighborhood of the naked singularity. Therefore the power-law 
dependence in the LTB space-time described above is basically realized by the analytical discussion about the case of 
the Minkowski space-time. Further the propagation effect due to the existence of the dust cloud and the occurrence 
of the central naked singularity is rather small. In other words, there is no peculiar behavior of ips even in the 
neighborhood of the central naked singularity. However, it should be noted that the odd-parity metric perturbation 
does not vanish in the neighborhood of the central naked singularity and Cauchy horizon associated with it. As a 
result, we conclude that the central naked singularity formation in the LTB space-time is 'marginally' stable against 
the odd-parity metric perturbations. 

We note that our analyses are not sufficient to determine the stability of the naked singularity formation in the 
LTB space-time. There remain some problems to complete the analysis. The first problem is to take account of 
odd-parity matter perturbations. We are now investigating this problem. The second is to consider the even-parity 
mode in which the metric and matter perturbations are essentially coupled with each other. This problem will be 
analyzed in future. In this paper, we have dealt with the marginally bound case. For the case that is not marginally 
bound collapse, the condition of the appearance of the central naked singularity is slightly different from the above 
case [ p3p^ | and hence there is a possibility that the behavior of ip s in this case is different from the marginally bound 
one. This case is now under investigation. 
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APPENDIX A: ANALYTIC SOLUTION OF THE DIPOLE MODE 



We present an analytic solution of the wave equations for the I — 1 mode. Substituting I — 1 into Eqs.( 3.16 ) and 
(3.17), we get 



ft(izV«) =o, 

d r [R^ips - 167rJ(r)] = 0. 



Eq.(A2) is easily integrated and we obtain 



- 16tt J(r) = c(t), 



where c(t) is an arbitrary function of t. Substituting this equation into ( Al), we obtain 

d t c{t) = 0. 



(Al) 
(A2) 

(A3) 

(A4) 



Hence the function, c(t), is temporally constant. From Eq.(3.20) and the regularity conditions for R, Lq and tp s , 
the left hand side of Eq.( |A3[ ) vanishes at the center. Therefore the function, c(t), should vanish identically and the 
solution for the dipole mode is obtained in the form 



*. = ^J(r). 



(A5) 



From the regularity condition, ( 3.2S| ), J(r)/R 4 , is identically zero at the regular center and hence this mode does not 
affect the formation of the central naked singularity. 
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FIG. 1. Conformal diagram of the LTB space-time with a globally naked singularity. ) denotes future (past) timelike 

infinity respectively, while i denotes spacelike infinity. £f + {£[~) denotes future (past) null infinity respectively. The dotted 
line H + indicates a future Cauchy horizon associated with the central naked singularity. The broken line is a null hypersurfacc 
on which we put initial wave packets. The initial positions of the wave packets are classified into Cases 1-3. For the locally 
naked singularity case, the Cases 1-3 are defined in the same manner as the globally naked case. 
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(a) Casel,2,3 (b) Case2 

FIG. 2. Plots of the gauge-invariant variable, ip 3 , at the center, r = 0, with an initial width a = 0.05rt for the Cases 1-3 for 
globally naked cases. In (a), the dotted line denotes Case 1 (r' c = 0.5rt), the broken line shows Case 2 (r' c = 0.82r&), and the 
solid line denotes Case 3 (r' c = 1.2r 6 ). In (b), the results of Case 2 are shown in more detail. The broken line in (b) is the same 
as the broken line in (a). The solid and dotted lines show the cases that the wave packets were put on the initial surface at 
r' c = 0.76rb and at r' c = 0.88n,, respectively. 
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FIG. 3. Plots of the gauge-invariant variable, ip s , at the center, r — 0, with an initial width a = 0.02r(, for Cases 1-3 for locally 
naked cases. The dotted, broken, and solid lines denote Case 1 (r' c = 0.28r b ), Case 2 (r' c = 0.38r 6 ), and Case 3 (r' c — 0.5Sr b ), 
respectively. 
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FIG. 4. Plots of the gauge-invariant variable, ip s , at the center, r = 0, in the LTB space-time for various widths of initial 
wave packets. The wave packets are put at r = 0.88rb on the initial null surface. The widths of initial wave packets are varied 
from 0.03r(, to 0.12r(,. The solid line (a) corresponds to the wave form of the wave packet with the initial width a — 0.03rv 
The broken line (b) is that of the initial width a = 0.05r;, while the dotted line (c) corresponds to that of the initial width 
a — 0.087V The broken dotted line (d) is that of the initial width a — 0.12r6- 
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FIG. 5. The relation between the widths of initial wave packets and the maximal values of \ip a \ at the center, r = 0. The 
results for the case of the LTB space-time with globally naked singularity are marked by open circles. The results of the locally 
naked case are marked by the triangles. The results of the Minkowski space-time are marked by cross marks. The broken line 
denotes the relation, \ip s \max oc ct~ 3 . 
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(c) a = 0.02 - 0.5r b in Case2 



FIG. 6. Plots of ips with an initial width a = 0.05rt at R = 100M as a function of time, t, in the LTB space-time. In (a), the 
solid line shows the result for the case with an initial time when we put the wave packet on the initial surface U/M = —65.310 
(Casel), the broken line is for U/M = -38.529 (Case2), the dotted line is for U/M = -13.610 (Case3). (b) depicts the details 
of the Case2. The solid line shows the plot of ip s with the initial time, ti/M — —34.336, the broken line is for U/M = —38.529, 
the dotted line is for U/M — —44.677. We find no diverging tendency of the gauge-invariant ip s when it approaches the Cauchy 
horizon. In (c), we vary the width of the initial wave packet in Case2. The solid line is a plot of a = 0.02rt, the broken line is 
the case of the initial width a — 0.25r6, the dotted line is that of a = 0.5rt. 
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(a) Case 1 (b) Case 2 

FIG. 7. Results of the comparison of the wave forms at the center. In (a), the solid line shows the LTB case, as a function 
of t/M, that is the same as the Case 1 in Figj^(a). The dotted line shows the corresponding one of the Minkowsi case, as a 
function of t, where a = 1.18 and r c = 13.9. In (b), the LTB case is the solid line and Case 2 in Fig.^(a). The Minkowski case 
is the dotted line where a = 1.25 and r c = 21.4. 
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FIG. 8. Wave forms along the constant circumferential radius both for the Minkowski and LTB cases are plotted. The 
solid line shows the LTB case, as a function of t/M, that is identical with Case 1 in Fig|] (a). The dotted line shows the 
corresponding one of the Minkowski case, as a function of t, where a = 1.33 and r c = 100. 
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